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We show that the spectrum of the Schrodinger equation in two or higher dimensions does not 
change when Dirichlet boundary conditions are enforced on a number of isolated points inside the 
original domain (piercings). We have obtained the analytical solution for spherically symmetric 
state of the d-dimensional simple harmonic oscillator pierced at the origin. Results for the case with 
multiple piercings are obtained numerically and agree with the theoretical prediction. In the case of 
a two dimensional parabolic quantum dot with two electrons and a single piercing in the origin we 
show that the energy of the ground state calculated to first order in perturbation theory goes over 
to the equivalent result in absence of piercing as the radius of the piercing becomes infinitesimal. 
Interestingly, we find that the leading finite size correction to the interaction energy is negative while 
the corresponding correction to the single particle energies is (as expected) positive. For a finite 
radius of the piercing and above a critical coupling, the first term dominates the second and the total 
energy of the dot with piercing is lower. Unfortunately the critical coupling found for this example is 
nonperturbative. Finally we study configurations allowing bound states in the continuum, showing 
that bound states survive to the insertion of piercings and that their energy is unchanged when the 
radius of the piercing vanishes. 

PACS numbers: 03.30.+p, 03.65.-w 



I. INTRODUCTION 



It is a well known result that the frequencies of a membrane do not change when Dirichlet boundary conditions 
are enforced on a number of isolated points inside the domain of the membrane: this curious result, which at first 
seems to be at odds with common sense, was first conjectured by Lord Rayleigh longtime ago More recently 
this problem has received attention in a number of papers, see for example 0, 0, I3j Hi @j 0il3l- In particular, in 
ref. Q Wang has obtained a general expression for the correction to the frequencies of a membrane with N internal 
circular regions of infinitesimal radius where no vibration occurrs: these corrections are seen to decay proportionally 
to l/logrg, ro being the radius of the circular cores. In the present paper we show that these interesting results can 
be extended to Quantum Mechanics, at least for non-interacting particles. In Section [TT] we derive a general formula 
for the correction to the energies of the Schrodinger equation in any number of dimensions, d > 2, with N internal 
points of infinitesimal radius, where Dirichlet boundary conditions are enforced (we will refer to these points as to 
"piercings" throught all the paper). The case d = 2 confirms the logarithmic behaviour of the corrections already 
observed in the classical membrane vibration countepart, while a power-like behaviour is found for higher dimensions. 
In Section Hm we use the simple harmonic oscillator with a single piercing in the origin (pSHO) as the natural testing 
ground of the general results obtained: in this case we are able to confirm analytically the predictions of Section 
im The case with multiple piercings, on the other hand, is studied numerically, leading to a further confirmation 
of the results. An interesting question that may be asked is whether the same conclusions of Sections [ll] and IIIII 
can be reached for N interacting quantum particles in d dimensions: since the problem is clearly too difficult to be 
attacked on general grounds, we have restricted our attention to a two dimensional parabolic quantum dot, with two 
electrons, which has been studied in ref. @. In this case we have been able to show that, to first order in perturbation 
theory, the energy of the 2D quantum dot with a piercing in the origin is unchanged for —^ 0; on the other hand 
we have obtained that the leading order finite size correction to the energy is not definite positive and indeed may 
become negative if the coupling of the Coulomb interaction is strong enough. Although the critical coupling where the 
change in sign occurrs is too large to trust a first order pcrturbative calculation, we believe that this result opens the 
question on whether the energy of the ground state of a quantum dot could be lowered adding one or more piercings 
of small but finite size. In Section |V] we have applied our results to the study of open systems supporting a bound 
state [13, [m [H, [H, , obtaining numerical results for the crossed wire configuration of Ref. [ll| which confirm our 
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predictions. In this way we prove that these bound states survive to the inclusion of piercings, even for piercings of 
finite extension. Finally, Section [VII contains our conclusions. 



II. GENERAL CONSIDERATIONS 



Consider the stationary Schrodinger equation (SSE) in d spatial dimensions 

^2 

- —Vl-ip{xi,. ..,Xd)+ V{xi, . . .,Xd)'4'ixi, ...,Xd)= Eipixi, ...,Xd) (1) 
im 

where {x,y) G V and Dirichlet boundary conditions are enforced on dV. Let Pi = {x^^\ . . . ,a;^'^) be a point internal 
to V, where Dirichlet b.c. are also imposed. The index i ranges from 1 to N, i.e. the total number of internal points 
with Dirichlet be. We call dV* = dV U J^f^i Pi- In this case the SSE reads 

- — Vd<^(a;i, ... ,Xd) + V{xi, . . . , Xd)4'ixi, . . .,Xd) ^ E'(j){xi, . . . ,Xd) , (2) 
2m 

E' and 4> being the new eigenvalue and wave function respectively. Assuming that V{xi, . . . , Xd) is real and using the 
two equations we can write 

{E - E'W^ = V^V- - VV^r ] - V • J , (3) 



where 



2m 



*Vd^-^Vd0* . (4) 



After integrating eq. ^ over I? one has 



E' -E ^'^^ ' '^'^'^'^ (5) 

/x,0*(2;,y)V'(2;,2/)rf''a; ■ 

We use Gauss theorem to convert the integral over the d-dimensional volume of D in the numerator into an integral 
over dV*: 

V-Jd'^x^ J-dSd-i + y^ J-dSf_^, (6) 
"D Jav .^^^ J Pi 

where Sd-i is the surface of T) and S^_-^ is the surface of a d-sphere centered in Pi with infinitesimal radius. Since 
both ij) and (j> obey Dirichlet be on dV, the first integral vanishes. 

On an d — 1 sphere of infinitesimal radius tq around the point Pi where Dirichlet boundary conditions are imposed 
one can express locally the solution, as a linear combination of the regular solution ^ and of a solution ^ which 
diverges at P,;: 

w V + WdC (7) 

where Vd is a constant to be determined by the condition 0(ro) = 0. 

Keeping only the laplacian term in the Schroedinger equation, we approximate ^ with the solution of the d- 
dimensional Poisson equation, diverging at Pi. 



We therefore obtain 



C V'('-o) d = 2 
Vdirn) = S '°s'-o (9) 



3 



After taking into account these observations, and neglecting the terms which vanish after integration we have 



V • JcTx 



■D 



2m 



E 



Pi 



N 2 



i=l 



(10) 



where S^_^ is the surface of a d — 1 sphere, which is given by 

2^d/2 d-l 

T{d/2) 



S, 



Therefore we may write: 



V • Jd^x 



(i) 

T{dj2) 



d = 2 
d>2 



(11) 



(12) 



Let us now come to the denominator in eq. ([5]): since we expect that the wave function in the presence of piercing be 
sensibly modified only in the close neighborhood of the points themselves, and the wave functions being normalized, 
we may argue that (p* {x, y)ip{x, y)d'^x ~ 1. The leading behaviour of the energy for tq ^ can thus be obtained 
from eq. 

An alternative -qualitative - argument to convince the reader of the isospectrality could be the following: consider 
the wave function 



,Xd) = ipi^i, ■ ■ ■ ,Xd) S(a;i, . 



,Xd) 



where 



N 



E{xi, ...,Xd) = ^'^^^''(^i' ---^^d) 



(13) 



(14) 



and is the solution to the homogeneous Poisson equation singular at a point Pi internal to D and vanishing on 
an infinitesimal circle of radius rp. ip{xi, . . . , Xd) is the solution in absence of piercing. We substitute cj) in eq. ^ and 
we obtain the equation: 



2m 



2^'iP ■ VS + i/^V^S 



+ (E-E')<j>{xi,...,xd)=0 



If we multiply this equation by (f) and integrate over the d dimensional volume we have 



iE-E') = 



2m 



2VV' • VS + i'V^E 



dS 



(15) 



(16) 



Since the expression in the square brackets is a function which is singular at all the internal points -P;, we expect 
that the main contribution to the integral comes from these regions. Around these points, however, = 0, and 
therefore we may expect that 



2VV' • V5 + V'V-^S 



d'^x « , 



(17) 



implying that E' E. 

We can also provide a more direct proof of this statement by considering the special case of a single piercing in two 
dimensions, which can then be generalized to N piercings in d dimensions. In this case eq. ()13p reads 



(x,y) = il^ix.y) log — 



As it stands this expression is not normalized and we should rather write 



\J lV'(a;,2;)plog^ 7^rfa;dy 



(18) 



(19) 
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If wc now take the limit tq — > cx3,for r > tq, we have 

hm 4>{x,y) = ip{x,y) (20) 

which supports the conclusions previously reached. 

Since the results obtained in this Section concern the inclusion of isolated Dirichlct points in the d dimensional 
domain where the SSE is defined, the reader may wonder at this point if it is possible to find regions of dimension 
< d' < d — 2 which leave the spectrum invariant when Dirichlet be are enforced. In this case the application of 
Gauss theorem must be specific to the problem, since >5'^_x will be now the infinitesimal surface enclosing the Dirichlet 
region. We may however convince ourselves that it is indeed possible by looking at a specific example: consider the 

Hamiltonian for a d > 3 problem, which is separable in the form H ~ + Hk+i d, where k > 2. Here i/i...fc 

contains the coordinates xi, . . . ,Xk and their conjugate momenta, while Hk+i....,d contains the remaining coordinates 
and momenta. In this case the solution of the problem is obtained as the direct product of the solutions of and 
Hk+i,...,d respectively. Therefore, adding N piercings of infinitesimal size to Hi ^ can be interpreted as adding N 
regions of dimension d — k with Dirichlet boundary conditions to the original domain. Since the spectrum of iJi...^ is 
unchanged under this operation, we conclude that the same is true for the total hamiltonian H. 



III. THE SIMPLE HARMONIC OSCILLATOR IN d DIMENSIONS PIERCED IN THE ORIGIN 



The simple harmonic oscillator in d > 2 dimensions is the ideal test of the results presented in the previous section, 
given that analytical results can be obtained. The SSE in this case reads 



2^ 



Vl + lfi^\xl + --- + xl)~E 



Tp(xi, ...,Xd) = , 



where is the Laplacian operator in d dimensions. For simplicity we study only spherically symmetric states with 

zero angular momentum and therefore consider only the radial part of the Laplacian (V^ + ^-jr^-§p)- ^ 

After introducing the standard dimensionless variables p = ^uj/Ti r and e = E /fiuj ^ we write "(/'(p) = e"'' ^^£,{p) 
and obtain the differential equation 



C'ip) 



2pU'(p) + (2e-d)e(p)=0. 



(21) 



We wish to solve this equation subject to the conditions ^(po) = and limp_,oo tpip) ~ 0, so that the wave function 
be square integrable on the domain. After enforcing the first condition we obtain the solution: 



ap) 



u 



d-2e d 



U 



d-2e d 



^c, [Po) 



(22) 



where U and L are the hypergeometric U and the Laguerre L functions. Since the Laguerre function blows up at 
infinity the full solution also blows up unless the condition 



U 



d-2e d 



= 



is met. As we see in Fig. [U which displays il{£, po) 



' 2" 

c/(^.i.pg) 



(23) 



, — ;= — for On = 10 ^, the zeroes of this function are 
logio i/po r(f ) 

slightly displaced to the right with respect to the unperturbed values, represented with circles. Notice that for specific 
values of the energy, the expression for the wave function simplify: for example, in two dimensions, for e = 2k + 1, 



^ Notice as well that the remaining states authomatically comply with the Dirichlct boundary condition at the origin and therefore both 

their energies and wave functions are unaffected. 
^ Actually, since the numerical results discussed in this paper are obtained for unit mass and sotting h = u) = 1, the use of p or r is 

equivalent. 
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with k integer, one goes over to the wave functions of the excited states of the unpierced osciUator (this corresponds 
to setting po in the largest positive node of such wave function). For e = 3 one has ^(p) oc (1 — p^), corresponding to 

Po = 1; for e = 5, one has ^(p) cx (2 — 4p^ + p^), corresponding to po = \/2 + \/2. 

In particular we are interested in knowing the behaviour of the energies of the ground state for po — > 0: for d = 2 
we obtain that e^^^ w 1 - 1/logpo, whereas for d > 2 we have e^'') w d/2 + c''') p^"^, where c^''' = 2/r(d/2 - 1). 
Notice that this expression for c'*^^ coincides with the general expression in eq. (|12p . 

We can now use eq. (1^^ to extract the dominant behaviour of the exact solution as po 0. In the limit r ^ 0, for 
the ground state, the confluent hypergeomctric function behaves as 

, d f — 2alogr d — 2 , ^ 

Uia,-,y)^[^rj^ d>2 (24) 

which confirms eq. (|13p . 

In the left plot of FigOwe display the quantity S(ro) = (ed — d/2)/rQ^'^ for the ground state of the SHO pierced 
in the origin in d = 2,3,4,5 dimensions for small values of rg. The horizontal lines are the coefficients obtained 
directly from the exact solution and from the general formula, given in the previous section. In the right plot of Fig[2] 
we compare the exact solution of eq. (|22p and the asymptotic one of eq. (|13p . both corresponding to a piercing radius 
ro = 10"^°, with the solution for the oscillator without piercing. 

The energy of the ground state of the pierced sho can be estimated precisely even in the case where the radius of 
the pierced region is sizeable using a variational approach. In this case we use a trial wave function 



for d = 2 and 



V'(p) = N e-p"'^ log (^^^ (1 + ap"") 



(25) 



(26) 



for d > 2. Here is a normalization constant and a and b are variational parameters. In the left part of Fig. [3] we 
show the energy of the ground state of the pierced sho, as a function of the pinning radius, tq, for d = 2, 3, 4, 5. In the 
limit ro ^ the energies tend to the energy of the unpierced sho, whereas for ro > 1 the curves approach a universal 
behaviour, which is dominated by the potential energy. Notice that the symbols correspond to the energy calculated 
numerically using the exact expression, while the lines correspond to the variational results. 

We briefly address the problem of a SHO with N piercings in two dimensions. Although this problem cannot be 
solved exactly for an arbitrary radius of the piercing, its spectrum should still go over to the unperturbed spectrum as 
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TABLE I: Theoretical and numerical results for the energy and the leading coefficient c'^-* of a SHO which contains A'^ piercings 
uniformly distributed on a circle of radius 1/3. 



N ejh epiT c^feViV c^j^/N 

2 1 0.9999977 0.8948766 0.8950938 

3 1 0.9999314 0.8948683 0.8987901 

4 1 0.9999960 0.8949873 0.8950401 
10 1 0.9998796 0.8956917 0.8976956 
20 1 0.9993188 0.8983882 0.9036415 



[I] 
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FIG. 2: Left plot:H(ro) = (ed — d/2)/rQ^^ for the ground state of the simple harmonic oscillator in d dimensions, pierced at 
the origin. Right plot:Wave function for the ground state of the sho: the solid and dashed lines correspond to the exact and 
asymptotic solutions for the sho pierced at the origin; the dotted line is the unpierced sho. The parameter ro = 10^^" is used 
in the first two curves, (color online) 




the radius of the piercing goes to zero, with a logarithmic strength which is given by eq. (|12p . In Table U we compare 
the numerical results for the energy and the leading coefficient c^^\ obtained by fitting the energy obtained using 
the expectation value of the hamiltonian in the asymptotic wave function for ro going from 10^"^° to 10^^", with the 
theoretical predictions obtained in Section |TT1 The overlap between the asymptotic and unperturbed wave functions 
has been calculated numerically for vq = 10^'"''': this overlap, as anticipated, is fairly close to 1 and simply assuming 
it 1 would lead to c^^) = N g-i/^ w 0.894839iV. 

The right plot of Fig. [3] displays the energy of the ground state of the SHO with piercings uniformly distributed 
on a circle of radius R = 1/3. The plot of the asymptotic wave function for the SHO with 20 piercing of radius 
ro = 10^^° is shown in FigHl 



IV. GROUND STATE OF A TWO DIMENSIONAL PARABOLIC QUANTUM DOT 

We now consider a two electron system confined by a 2D parabolic potential in a zero magnetic field. Ref. contains 
an analytical expression for the first order contribution to the ground state energy of this system in nondegencratc 
perturbation theory. We will first review the main steps of this calculation for the unpierced quantum dot and then 
extend the calculation to the quantum dot pierced in the origin. 

The hamiltonian for this problem is 

^9 2 ? "2 9 9 9 



f f ' ' « f t » « 



_LLLL 



FIG. 3; Left plot: Energy of the ground state of the d-dimensional simple harmonic oscillator as a function of the radius ro of 
the internal region where the wave function vanishes. The symbols correspond to the numerical solution of the exact equation, 
whereas the dashed lines correspond to the variational estimate. Right plot:Energy of the ground state of the two dimensional 
quantum dot with A*' piercings uniformly distributed at a distance d — 1/3 from the origin, calculated using the asymptotic 
wave function. We use h — m — uj = 1 and plot the energy as a function of the radius ro. (color online) 





FIG. 4: Asymptotic wave function for the SHO pierced on 20 points uniformly distributed on a circle of radius R = 1/3. In 
this case ro = 10^ (color online) 
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where the last term provides the Coulomb repulsion between the electrons. 

In the absence of this term the total wave function is obtained as the direct product of the single particle wave 
functions of a simple harmonic oscillator, which are given by: 



/2n 



(28) 



where i™ are the associated Laguerre polynomials. Nnm = \j2n\cP- l(n + \m\)\ is the normalization factor and 

(29) 



a = \J [iLojU. The single particle energies are given by 

Enra = hLj{272 + |m| + 1), 

with n = 0, 1, . . . and m = 0, ±1, ±2, .... The single particle wave function for the ground state is simply 

$oo(r,0)-^e-"''^'/^ (30) 



In the ground state of the two particle system the electrons form a spin singlet and therefore the orbital wave 
function is symmetric in the electron coordinates: 



^'(fi,?^2) = $00 (^1)^00(^2) 



X+.1X-.2 - X+.2X-,1 
V2 



(31) 



Treating the Coulomb repulsion between the electrons as a perturbation, the authors of Ref. have obtained an 
analytic expression for the first order correction in nondcgencrate perturbation theory: 



J2„ / j2„_„-a^(r?+r^) ^ 



d ri d r2e 



— — — e a\i ^ , 



using the identity 



In - r2 



+00 „oo 



(32) 



(33) 



We may now discuss the same problem in the presence of a piercing in the origin. As previously found, the single 
particle wave function for the ground state in presence of piercing of infinitesimal radius rg is simply given by 



$oo(?^) = ^00 e-"''-'/2 log 



ro 



The normalization constant is expressed in term of the Meijer G function as 



\ 



IT G^;° a2ro2 



1,1 
0,0,0 



(34) 



(35) 



and N, 



00 



07 



for 7'o ^ 0. 



yTTlogaro ^ 2^7? log^ aro 

The perturbative correction to the ground state energy is therefore given by 



A^(i) = e'N^^ [ d\i [ d\2 log^^log^^ e-"'('-i+'-') 
J J ro ro 



e^A^o^o / dk 



rf2^1og2_e-" ^ jQ{kr) 



(36) 



^ Because the ground state is spherically symmetric only the m = term in this expression contributes. 
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In order to evaluate the integral in the square bracket we use the series representation 



3=0 



and evaluate the integral 



POO 

Ij = dr r^J+i log^ 

J rn 



(37) 



1,1,1 
0,0,0,j + l 



log2(aro) - log(aro)V'^°n.? + 1) + 



where + 1) is the polygamma function given by ip^^^z) = X^fcLi 

together with eq. ([57]) . to write eq. ([55)1 as 



fe fe+2-l 



(38) 

7. We use this expression, 



log'^aro-2 r 0, 



log(-^) ) log'Ko) 



TT K — \/27r7 
2 log aro 



(39) 



where k = -2 /„ dfce"— [T [Q, -xj + log ("xj j ~ 2.240700135. This equation tells us that in the limit jq ^ 
the interaction energy to first order in PT is exactly the same as in the case of the quantum dot without piercing, 
while the leading finite size correction to this energy is negative. In other words, at least to first order in PT, the 
presence of a piercing lowers the interaction energy. Combining this expression with the leading order expression for 
the single particle energy we obtain the total energy to first order in PT goes as 



Etot ^ huj {2 + X\/^ - 

2 log aro 



2 - (k - V27r7)A 



(40) 



where A = ^-j , using the notation of Ref . Q . For A > A,; 



= w 2.51941 the total energy of the system is 

lower in presence of a piercing, for sufficiently small values of tq. On the other hand we know that for large tq, the 
single particle energy grows rapidly, thus suggesting the presence of a minimal energy at a finite tq. Fig. [5] confirms 
this expectation and shows the energy of the ground state of the two dimensional quantum dot with piercing for A = 3 
calculated to first order in perturbation theory, using the single particle wave functions of eq. (|22p . The interaction 
energy in this case is evaluated numerically. 

Notice that a fit of the numerical results obtained for 10^^" < ro < 10^^" yields 



4'/t(^o) 



1.25322 



0.79415 0.210805 0.572547 



log aro log aro log aro 
which agrees remarkably well with the theoretical result previously obtained 



Eg\ro) = e'a 



27r7 



log aro 



e^a 



1.25331 



0.793835 
log aro 



(41) 



(42) 



V. BOUND STATES IN THE CONTINUUM 



We wish now to apply the general results obtained in Section |TT] to the case of the Helmholtz equation on domains 
which extend to infinity, but consist of wires with crossings and/or bendings. It is well known that in this case the 
spectrum of the Laplacian contains one or more bound states, depending on the number of crossings and bendings, 
below the threshold of the continuum. For example, Exner and Seba first showed in 10| that a smoothly curved 
waveguide can support a bound state; Schult et al. have reached a similar conclusion in [ll| for a different configuration 
consisting of two crossed wires, of infinite length. Avishai and collaborators have also proved the existence of a bound 
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10""° 10"" 10"'' 10"' 10""* 10° 



FIG. 5: Energy of the ground state of the two dimensional quantum dot with pinning to first order in perturbation theory for 
A = 3 and 7i = ?n = a; = l, asa function of the radius ro. The horizontal line is the equivalent result in the absence of pinning, 
(color online) 



state in the broken strip configuration for arbitrarily small angles (see [l^), while Goldstone and JaflPe [13[ have given 
a variational proof of the existence of a bound state for an infinite tube with bendings in two and three dimensions. 

The general results obtained in our Section |TT] suggest that when a number of piercings of infinitesimal radius is 
added to a crossed wire or to a bent waveguide, the energy of the bound state remains precisely the same although 
the corresponding wave function changes. If the radius rp of the piercing is now made finite, still one expects the 
bound state to survive, up to some critical value of tq. Notice that the considerations made at the end of Section HIl 
also indicate the possibility to extend the same results to filaments inside three dimensional bent tubes. 

In this Section we provide an explicit confirmation of our general results, by studying the crossed wire configuration 
of Ref . [m . In analogy with Ref. [ll| we have decided to study this problem using two methods: the first method 
is based on a collocation approach which allows one to obtain a numerical solution to the Helmholtz equation; the 
second method uses an expansion in a complete set of solutions of the Helmholtz equation on the five domains which 
compose the global domain (shown in Fig.l of Ref. [Tll| ) . The explicit expressions for these sets of functions may be 
found in the paper of Schult ct al. 

Let us first briefly describe the first method. In this case we have discretizcd the Helmholtz equation on a uniform 
grid using a set of functions, called Little Sine Functions (LSF), derived in Ref. [l^. These functions have recently 
been applied in [l6j to the study of vibration of membranes of arbitrary shapes. Although we refer the reader to 
Ref. [l6| for the technical details concerning the implementation of this method, we just mention that the inclusion 
of piercings inside the domain is handled with extreme simplicity in this approach, provided that the piercing falls on 
one of the point forming the mesh. As a matter of fact, each LSF is an approximate representation of a Dirac delta 
function, peaked on one of the points of the mesh: the exclusion of a point from the mesh is therefore obtained by 
eliminating the corresponding LSF from the set of functions used to discretize the Hamiltonian. 

Using the collocation method we have solved numerically the Helmholtz equation on the crossed wire domain, 
assuming that the arms have a width w = 2 (we have also set ?i = m = 1 in our calculation). Although the arms 
have infinite extension, we have used arms of finite length L = 10, in other words we have limited the crossed wire to 
the interior of a square of sides 20 x 20. This choice is expected to affect very mildly the lowest energy state which is 
localized at the crossing. The energy obtained in this way provides in any case an upper bound to the true energy. 

In Fig. [S] we display the energy of the ground state obtained numerically with the LSF, for the case of a crossed wire 
without piercings (circles) and for a crossed wire with a piercing in the origin. The horizontal line is the threshold 
energy for the continuum. These results are obtained for specific values of the grid spacing, h, for which the boundary 
of the crossed wire falls precisely on the mesh: as discussed in Ref. [iGj , in this case one obtains a monotonous sequence 
of values which converge to the exact result from above. We have extracted the continuum limit fitting the two sets 
of data with 



(43) 
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FIG. 6: Energy of the ground state of the crossed wired of width w = 1 calculated using the collocation method as a function 
of the grid spacing h. The circles are the results for the crossed wire without piercing and the crosses are the results for the 
crossed wire pierced in the origin. The horizontal line is the threshold for the continuum, (color online) 



TABLE II: Ground state energy of the crossed wire obtained with the collocation method. Eq jg ^jjg energy of the 

configuration without (with) piercing. The last row are the energies obtained using eqs. (|43|) and (|44|) with a least square 
approach. 
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for the case without piercing and 



(,) _ l+h/l0gh + b2h>'^ 

- ^" l + b./logh + hhOo ' ^^^^ 



for the case with piercing. The coefficients and bi are obtained by a least square procedure. The last row of the Table 
contains the results obtained from the coefficients ao and bo of the fits above using a least square procedure: the results 
for the two cases are almost degenerate and agree with the result quoted in Ref. [ll|, i.e. kw = "^^'"'^ ~ 0.8127r. 
Let us now describe the second method. As done in [ll| we express the solution in each of the five regions as a 
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FIG. 7: Energy of the ground state of a crossed wire with a piercing in the origin as a function of the radius ro. The horizontal 
line is the threshold of the continuum, (color online) 



linear combination of elementary solutions fulfilling the Helmholtz equation: 
M^,y) = v\„cosh^e-(-/-^) COS 

n 2 W 



N-1 



(2n+l)ny , (2n+l)7ra; 

cosh PnX COS h cosh PnU COS 



w 



(45) 
(46) 



where pn = y '''^"^^■i ^ ^"f'" . Notice that there is no need of writing the remaining solutions, since they are 

obtained from il)i{x,y) by means of rotations and reflections. The unknown coefficients 6„ are obtained by imposing 
the continuity of the normal derivative of i}}[x,y) at the border between the regions I and V. Using iV = 8 we have 
obtained E = 0.849555 (using as before w = 2 and ?i = ?n = 1), corresponding to kw = « 0.8298337r, which 

is about 2% larger than the result previously obtained (this is consistent with the result found in (l^). 

As we have seen previously the inclusion of a circular piercing is obtained by multiplying the original solution by 
a factor log r/rg, rg being the radius of the piercing. If this is correct, one should see that, as is sent to zero, the 
energy approaches the value in the absence of piercing. 

In Fig. [7] we show the energy of the ground state of the crossed wire with a piercing of radius ro in the origin as a 
function of vq itself. Using once more a least square procedure we find that in the limit rp — > 0, £'o ~ 0.848, which 
is remarkably close to the value without piercing. The critical value of the piercing radius for which the threshold 
energy is reached is rg""'* « 0.0899. In Fig. [5] we show the wave function of the ground state of the crossed wire with 
a piercing of radius rO ~ 10^^ at the origin, obtained using the expansion in terms of the elementary solutions. 



VI. CONCLUSIONS 



In this paper we have proved that the spectrum of the d-dimcnsional Schrodinger equation does not change when 
piercings of infinitesimal size are added to the d-dimensional domain. The example of the simple harmonic oscillator 
in d dimensions is worked out and the expected results are obtained both analytically (for a single piercing in the 
origin) and numerically (up to 20 piercings). Using these results we have considered a two dimensional parabolic 
quantum dot and we have calculated the energy of the ground state to first order in perturbation theory, up to the 
leading finite size correction in the piercing radius. We have found that the interaction energy of the quantum dot 
is lower for piercings of finite size, and that it can dominate the corresponding finite size correction to the single 
particle energy above a critical coupling. In our calculations the critical coupling turns out too be large to trust a 
first order perturbative result. This outcome should motivate a further study of this system, either involving higher 
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FIG. 8: Asymptotic wave function of the ground state of a crossed wire with a piercing of radius ro = 10 ^. (color online) 



order pcrturbativc corrections or a variational calculation, which we hope to carry out soon. Another application 
considered in the present paper is to configurations supporting bounds states in the continuum, such as wires with 
crossings and bendings. We have explicitly shown that the inclusion of piercings to these systems does not alter the 
energy of the ground state, contrary to naive expectations and in perfect accord with our general considerations. 

A further remark that we wish to make concerns the Casimir effect on domains with piercings: at least in the case 
of non interacting fields, and for piercings of infinitesimal size, our results imply that no net effect should appear when 
piercings are added to a domain. In presence of interactions, whose study is certainly a formidable task, the inclusion 
of piercings may affect the spectrum thus leading to a net effect. 
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